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We simulate the central reactions of nearly symmetric, and asymmetric systems, for the 
energies at which the maximum production of IMFs occurs (13^'^^'^). This study is carried 
out by using hard EOS along with cugnon cross section and employing MSTB method for 
clusterization. We study the various properties of fragments. The stability of fragments 
is checked through persistence coefficient, gain term and binding energy. The information 
about the thermalization and stopping in heavy-ion collisions is obtained via relative 
momentum, anisotropy ratio, and rapidity distribution. We find that for a complete 
stopping of incoming nuclei very heavy systems are required. The mass dependence of 



various quantities (such as average and maximum central density, collision dynamics as 
well as the time zone for hot and dense nuclear matter) is also presented. In all cases 
(i.e., average and maximum central density, collision dynamics as well as the time zone 
for hot and dense nuclear matter) a power law dependence is obtained. 

1 Introduction 

Heavy-ion collisions at intermediate energies offer an opportunity of producing nuclear 
systems under the extreme conditions of temperature and density. At high excitation 
energies, the colliding nuclei compress each other as well as heat the matter [1-3]. This 
leads to the destruction of initial correlations, as a result the matter becomes homogeneous 
and one can have global stopping. The global stopping is defined as the randomization 
of one-body momentum space or memory loss of the incoming momentum. The degree of 
stopping, however, may vary drastically with incident energies, mass of colliding nuclei and 
coUiding geometry. The degree of stopping has also been linked with the thermalization 
(equilibrium) in heavy-ion collisions. More the initial memory of nucleons is lost, better 
it is stopped. 

The fragmentation of colliding nuclei into several pieces of different sizes is a complex 
phenomenon. This may be due to the interplay of correlations and fluctuations emerging 
in a collision. Several studies, in literature, have been made to check the fragmentation 
pattern. Fragmentation pattern has been found to depend on the size of the colliding 
nuclei, incident energy as well as impact parameter [1,4-6]. Dhawan et al. [7] studied 
the degree of stopping reached in intermediate energy heavy-ion collisions. They found 
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that degree of stopping decreases with increase in impact parameter as well as at very 
high energies. They suggested that the light charged particles (LCPs) (2<A<4), can 
be used as a barometer for studying the stopping in heavy-ion collisions. As lighter 
fragments mostly originate from midrapidity region whereas intermediate mass fragments 
(IMFs) originate from surface of colliding nuclei and can be viewed as remnants of the 
spectator matter. On the other hand, Sood and Puri [8] studied the thermalization 
achieved in heavy-ion collisions in terms of participant-spectator matter. They found that 
participant-spectator matter depends crucially on the collision dynamics as well as history 
of the nucleons and important changes in the momentum space occur due to the binary 
nucleon-nucleon collisions experienced during the high dense phase. The collisions push 
the colliding nucleons into midrapidity region responsible for the formation of participant 
matter. This ultimately leads to thermalization in heavy-ion collisions. Vermani et al. [8] 
used rapidity distribution of nucleons to characterize the stopping and thermalization of 
the nuclear matter. They found that nearly full stopping is achieved in heavier systems 
hke ^^^Au-|-^^^Au whereas in hghter systems a larger fraction of particles is concentrated 
near target and projectile rapidities, resulting in a broad Gaussian shape. The lighter 
systems, therefore, exhibit larger transparency effect, i.e., less stopping. Puri et al. [9] 
studied the non-equilibrium effects and thermal properties of heavy-ion collisions. They 
found that the heavier masses are found to be equilibrated more than the lighter systems. 

Recently, Sisan et al. [10] studied the emission of IMFs from central coUisions of nearly 
symmetric systems using a 47r-array set up, where they found that the multiplicity of IMFs 
shows a rise and fall with increase in the beam energy. They observed that (the 
energy at which the maximum production of IMFs occurs) increases linearly with the 
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system mass, whereas a power-law (oc A"^) dependence was reported for peak multiplicity 
of IMFs with power factor r = 0.7. Though percolation calculations reported in that paper 
failed to explain the data, subsequent calculations using QMD model [11] successfully 
reproduced the data over entire mass. One is, therefore, interested to understands how 
nuclear dynamics behaves at this peak energy i.e. whether linear increase reported for 
the multiplicity of fragments remains valid for other observables or not. We here plan 
to investigate the degree of stopping reached and other related phenomena in heavy- 
ion reactions at peak center-of-mass energies. We also check system size dependence of 
average and maximum central density, collision dynamics as well as the time zone for hot 
and dense nuclear matter at peak center-of-mass energy. 

This study is made within the framework of the quantum molecular dynamics model, 
which is described in detail in Refs. [1,4, 11-18]. 

2 Results and Discussion 

For the present study, we simulate the central reactions (b =0.0 fm) of ^°Ne-|-^°Ne, 
^°Ar+45Sc, ^^Ni+s^Ni, 86Kr+93j^b, ^^^Xe+^^^Sn, ^^Kr+^^^Au and ^^^Au+^^^Au at the in- 
cident energies at which the maximal production of intermediate mass fragments (IMFs) 
occurs. These read approximately 24, 46, 69, 77, 96, 124, and 104 AMeV, respectively 
for the above mentioned systems [11]. Note that in lower energy region phenomena like 
fusion, fission and cluster decay are dominant [19]. Here we use hard (labeled as Hard) 
equation of state along with energy dependent Cugnon cross section {cr(j^'^) [15]. The 
reactions are followed till 200 fm/c. The phase-space is clusterized using minimum span- 
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ning tree method with binding energy check (MSTB). The MSTB method is an improved 
version of normal MST method [16]. Firstly, the simulated phase-space is analyzed with 
MST method and pre-clusters are sorted out. Each of the pre-clusters is then subjected 
to binding energy check [11, 16]: 

i=l ' j^i 

We take Ebind — -4.0 MeV if > 3 and -Bbmd = 0.0 otherwise. Here is the number 
of nucleons in a fragment and is center-of-mass momentum of the fragment. This is 
known as Minimum Spanning Tree method with Binding energy check (MSTB) [11, 16]. 
Note that nucleons belong to a fragment if inequality (1) is satisfied. The fragments formed 
with the MSTB method are more rehable and stable at early stages of the reactions. 

One of the important aspects in fragmentation is the stability of fragments as well 
as surrounding nucleons of a fragment. The change in the nucleon content of fragments 
between two successive time steps can be quantified with the help of persistence coefficient 
[17, 18]. 

Let the number of pairs of nucleons in cluster C at time t is xc{t)=0.5*Nc{Nc-l). At 
the time At later, it is possible that some of the nucleons belonging to cluster C have left 
the cluster and are the part of another cluster or are set free or others may have entered 
the cluster. Now, let Ncjy be the number of nucleons that have been in the cluster C at 
time t and are at t + At in cluster D. We define 

<^c{t + At)^J2 * ^CniNcn - 1)- (2) 
D 

Here the sum runs over all fragments D present at time t + At. 
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Figure 1: The persistence coefficient as a function of time for LCPs, HMFs, and IMFs. 
Dashed, dotted, dash-dotted and sohd lines are for ^"Ar+'^^Sc, '^^Kr+^^Nb, ^'^Kr+^^'^Au 
and ^^''Au+^'^^^Au, respectively. 



The persistence coefficient of cluster C can be defined as [17, 18] 



At. 



Pc{t + —) = ^c{t + At)/xc{t). 



(3) 



The persistence coefficient averaged over an ensemble of fragments is defined as: 



c 



At. 



(4) 



where Nj^ is the number of fragments present at time t in a single simulation. The quantity 



is then averaged over a large number of QMD simulations. The stability of a fragment 
between two consecutive time steps can be measured through persistence coefficient. If 
fragment does not emit a nucleon between two time steps, the persistence coefficient is 
one. On the other hand, if fragment disintegrates completely, the persistence coefficient 
will be zero. If we remove one nucleon from a fragment C, the persistence coefficient is 
Pc(t + At/2) = (Nc-2)/Nc i.e., 0.333 for Nc = 3 and 0.8 for Nc = 10. For example for 
mass 10, when one nucleon is emitted we have two entities at later time step consisting 
of free nucleon and fragment with mass 9. The Pc{t + Ai/2) is the contribution from 
all such entities existing at later times. It, then, measures the tendency of the members 
of given cluster to remain together. In fig. 1, we display the persistence coefficient 
for various fragments i.e., light charged particles (LCPs) (2<A<4), intermediate mass 
fragments (IMFs) (5<A<44) as well as heavy mass fragments (HMFs) (10<A<44). The 
various lines have been defined in the caption. It is clear from fig. that the saturation 
value of persistence coefficient is slightly higher in case of LCPs as compared to heavier 
fragments. One can conclude that the final fragments are formed after 130 fm/c when 
this coefficient is 0.8. Before that time there is a strong exchange of nucleons between 
the fragments. The number of medium and intermediate size fragments increases because 
the largest fragment falls finally into this mass bracket. The persistent coefficient reaches 
its asymptotic value later due to the interaction between fragments as well as between 
fragments and free nucleons. Due to this interaction nucleons are sometimes absorbed or 
emitted from the fragments. This process changes the details but not the general structure 
of the fragmentation pattern. 

The persistence coefficient tells about the stability of different fragments between two 
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Figure 2: The Gain term as a function of time for LCPs, HMFs, and IMFs. Lines have 
same meaning as in fig. 1. 

successive time steps. But it does not provide any information whether a fragment has 
swallowed some nucleons or not. To check this, we use a quantity called "Gain" [18]. The 
Gain represents the percentage of nucleons that a fragment has swallowed between two 
consecutive time steps. Let be the number of nucleons belong to a fragment a at time 
t. Let N;^^ be the number of nucleons which were in cluster a at time t and are in cluster 
(3 at time t + At. The Gain is defined as: 



a 



Gaznit + At/2) = x ^^^^^^^ ^"^^^ (5) 
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Figure 3: The time evolution of binding energy per nucleon for LCPs, HMFs, and IMFs. 
Lines have same meaning as in fig. 1. 

7] = 0.0, 0.5, and 1.0 if N^^ < 0.5 N^, N^^ = 0.5 and N^^ > 0.5 N^, respectively. 
Naturally, a true Gain for a fragment a is only if its nucleons constitute at least half of 
the mass of new fragment f3. The Gain term will tell us whether the interactions among 
fragments have ceased to exist or not. 

In fig. 2, we display the gain term for LCPs, HMFs, and IMFs. As discussed earlier, 
the value of persistence coefficient is slightly higher in case of LCPs. Therefore, gain term 
will be smaller for LCPs as shown in fig. 2. As evident from fig., in case of heavier systems 
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the gain term has higher value because of the large nucleon-nucleon interactions. 

To check the stability of fragments, we display in fig. 3, the binding energy per nucleon 
as a function of time for LCPs, HMFs and IMFs. We find that even at 200 fm/c, small 
fragments are still not cold and take very long time to cool down. Whereas the heavy 
fragments are properly bound having binding energy per nucleon around -5 to -7 MeV. 

The quantities which are closely related to the degree of thermalization are relative 
momentum {Kr) and anisotropy ratio {Ra)- The average relative momentum of two 
colliding fermi spheres is defined as [3,9,20]: 

{Kn)^{\Pp{v,t)-PT{v,t)\/h), (6) 

where 

Here and pj are the momentum and density experienced by j*'* particle and k stands 
for either target or projectile and r refers to a space point in central sphere of 2 fm radius 
to which all calculations are made. The {Kr) is an indicator of local equilibrium because 
it depends on the local position r . 

The second quantity is anisotropy ratio which is defined as [3, 7, 9, 20]: 

(Ra) = ^ . (8) 

The anisotropy ratio {Ra) is an indicator of global equilibrium of the system because 
it represents the equilibrium of the whole system and does not depend upon the local 
positions. The full global equilibrium averaged over large number of events will correspond 

to {Ra) - 1. 
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Figure 4: The time evolution of (a) relative momentum and (b) anisotropy ratio. Lines 
have same meaning as in fig. 1. 

In figs. 4a and 4b, we display, respectively, (-ft'i?) and (i?a) ratio as a function of time 
for different system masses. The initial value of relative momentum increases whereas of 
the anisotropy ratio decreases with mass of the system since '^^^m increases with increase 
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in the system mass. It is interesting to see that the relative momentum is large at the start 
of the reaction, and finally at the end of the reaction, the value of {Kr) is nearly zero. This 
means that at the end of the reaction, the local equilibrium is nearly reached. However, 
the saturation time is nearly the same throughout the mass range. It is clear from the 
fig. 4b anisotropy ratio changes to a greater extent during the high density phase. Once 
the high density phase is over, no more changes occur in thermalization. Interestingly, 
the heavier nuclei are able to equilibrate more than the lighter nuclei. This is because of 
the fact that the number of collisions per nucleon for the ^^^Au+^^^Au reaction is larger 
than for the ^°Ar+^^Sc reaction. 

The rapidity distribution is also assumed to give information about the degree of 
thermalization achieved in heavy-ion reactions. The rapidity distribution of i*^ particle is 
defined as [7]: 

Here E(i) and p^i'') respectively, the total energy and longitudinal momentum of i*'* 
particle. Naturally, for a complete equilibrium a single Gaussian shape peak is expected. 
In fig. 5, we display the rapidity distribution of free-nucleons, LCPs, HMFs as well as 
IMFs. Rapidity distribution of all types of fragments indicate that heavier systems are 
better thermahzed as compared to hghter ones. For hghter nuclei, we get relatively flat 
distribution. It is quite clear from the figs. 5(c) and 5(d) that very light nuclei (^°Ar-|-*^^Sc) 
exhibit a peak at target /projectile rapidity indicating a non-equilibrium situation. Hence, 
for a complete stopping of incoming nuclei very heavy systems are required. For HMFs 
and IMFs, in case of ®^Kr-|-^^^Au, peak shifts towards left due to asymmetric reaction. The 
effect is more pronounced for different kinds of fragments as compared to free-nucleons. 
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Figure 5: The rapidity distribution, dN/dY, as a function of reduced rapidity, Y cm. I-beam- 
Lines have same meaning as in fig. 1. 

In fig. 6, we display the participant, spectator, and super-participant matter (obtained 
at 200 fm/c and defined in terms of nucleon-nucleon collisions) as a function of the total 
mass of the system. All nucleons that have experienced at least one collision are termed as 
participant matter. The remaining matter is counted as spectator matter. The nucleons 
having experienced more than one collision are termed as super-participant matter. As 
evident from the fig. 6, participant matter shows a nearly mass independent behavior ( r = 
0.15 0.016). The spectator and super-participant matter exhibit a power law dependence 
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Figure 6: The participant, spectator, and super-participant matter (defined in terms of 
nucleon-nucleon collisions) as a function of composite mass of the system. The dashed, 
dotted, and solid lines represent the the fits with power law A"^. 

with r = -0.46 0.039 for spectator matter and r = 0.46 0.06 for super-participant matter. 

In fig. 7, we display the system size dependence of the maximal value of average 
density {p""^^) (solid circles) and maximum density l^p^'^^^ (solid stars). Lines represent 
the power law fitting (oc A'^). The maximal values of {p""""^) and l^p^'^^^ follow a power 
law (oc A^) with r being 0.08±0.02 for the average density (p"^9) and 0.034±0.008 for 
maximum density (p™"^^) i.e., a slight increase in density occurs with increase in the size 
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Figure 7: The maximal value of the average density {p°'^^)max (upper part) and maximum 
density {p"^"'^)max (lower part) as a function of the composite mass of the system. The 
solid lines are the fits to the calculated results using obtained with minimization. 



The average is done over all space points on a sphere of 2 fm radius at center-of-mass. 
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Figure 8: The time zone for pavg > Po (upper part) and forpmax > Po (lower part) as a 
function of composite mass of the system. The sohd hues represent the the fits with 
power law A"^. 
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Figure 9: The total number of allowed collisions versus composite mass of the system. 
The solid line represents the fits with power law A'^. 

of the system. This is because, '^^cm increases with the size of the system. 

Apart from the maximal quantities, another interesting quantity is the dense zone at 
the peak energy. In fig. 8, we display the time interval for which pavg > Po (solid circles) 
and pavg > Po/2 (open stars). Again both quantities follow a power law behavior with r = 
0.15±0.05 and r= -0.04±0.06, respectively, for pavg > Po and pavg > Po/2. This indicates 
that the time duration for which p^vg is greater than the normal nuclear matter density 
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increases with the mass of the system. 

The system size dependence of the (allowed) nucleon-nucleon collisions (solid squares) 
is displayed in fig. 9. The results are displayed at 200 fm/c where the matter is di- 
luted and well separated. The nucleon-nucleon collisions increase with the system size. 
This enhancement can be parametrized with a power law proportional to A"^ with r = 
1.28±0.054. At fixed incident energy nucleon-nucleon collisions should scale as A. This 
has been tested by Sood and Puri [21]. Here power factor is greater than one since with 
increase in mass of the system E^.^^^ also increases. 

3 Summary 

In the present study, we have simulated the central reactions of nearly symmetric, and 
asymmetric systems, for the energies at which the maximum production of IMFs occurs 
(E^^^'^), using QMD model. This study is carried out by using hard EOS along with 
cugnon cross section and employing MSTB method for clusterization. We have studied the 
various properties of fragments. The stability of fragments is checked through persistence 
coefficient, gain term and binding energy. We observed that at 200 fm/c, small fragments 
are still not cold and they take very long time to cool down whereas the heavy fragments 
are properly bound. The information about the thermalization and stopping in heavy-ion 
collisions is obtained via relative momentum, anisotropy ratio, and rapidity distribution. 
We found that for a complete stopping of incoming nuclei very heavy systems are required. 
The mass dependence of various quantities (such as average and maximum central density, 
collision dynamics as well as the time zone for hot and dense nuclear matter) is also 
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presented. In all cases (i.e., average and maximum central density, collision dynamics 
as well as the time zone for hot and dense nuclear matter) a power law dependence is 
obtained. 
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